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Abstract. Additive symmetric Levy noise can induce directed transport of 
overdamped particles in a static asymmetric potential. We study, numerically and 
analytically, the effect of an additional dichotomous random flashing in such Levy 
ratchet system. For this purpose we analyze and solve the corresponding fractional 
Fokker-Planck equations and we check the results with Langevin simulations. We 
study the behavior of the current as function of the stability index of the Levy noise, 
the noise intensity and the flashing parameters. We find that flashing allows both 
to enhance and diminish in a broad range the static Levy ratchet current, depending 
on the frequencies and asymmetry of the multiplicative dichotomous noise, and on 
the additive Levy noise parameters. Our results thus extend those for dichotomous 
flashing ratchets with Gaussian noise to the case of broadly distributed noises. 
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1. Introduction 

The study of noise induced transport in anisotropic spatially periodic systems is a 
relevant subject for many problems in physics [U[2] and biology [3], and is also acquiring 
an increasing interest due to possible applications in the development of technological 
devices [2]. 

In general, the basic models consider a particle in an external spatially periodic 
potential with broken symmetry (the ratchet-potential), and subject to an additional 
signal or source of fluctuations which generates the nonequilibrium condition necessary 
for the emergence of directional transport [2]. Thermal noise is also usually considered, 
although depending on the nature of the nonequilibrium forces it may be not a 
fundamental ingredient. Ratchet models can be classified in two main classes according 
to the way in which the non equilibrium forcing affects the particle dynamics. Simple 
additive nonequilibrium forces lead to the so called rocking ratchets, while, when the 
signal enters as a multiplicative modulation of the ratchet potential, we speak of flashing 
ratchets. Systems combining the two kind of forcings can also been considered [I]. 

Recently, two simultaneous papers [SI E] have shown that a minimal setup for 
producing directional transport is obtained by considering a simple static ratchet 
potential and an additive white symmetric Levy [H|6l[7] noise as the only two ingredients. 
Here we will refer to such system as the static or non flashing Levy ratchet. The preferred 
direction of motion for such system is found to be towards the steepest slope of the 
potential. The effect was explained in [5] as a consequence of the large Levy jumps, 
which lead the particles to the flatter zones of the potential with larger probability than 
to the steeper zones. In the limit case in which the stability parameter a [21 [6] defining 
the Levy noise is set to be equal to 2, the Levy distribution becomes Gaussian[5j E] and 
the equilibrium situation with vanishing current is recovered. 

After the mentioned pioneering papers on the static Levy ratchet, several works 
have appeared providing further analysis of the system and studying different gener- 
alizations. In [8] it has been shown that an inversion of current can be obtained by 
considering a time periodic modulation of the chirality of the Levy noise. Reference [9] 
studies inertial effects and propose a way of measuring the rectification efficiency in Levy 
ratchets. In [10] . the weak noise limit is analyzed. In [UJ, the related problem of the 
spatially tempered fractional Fokker-Planck equation is studied. Reference [12] analyzes 
the competition between Levy forcing and a periodic a.c. driving in a ratchet system, 
while [13J studies the coexistence of Levy flights and subdiffusion. The increasing inter- 
est in ratchet systems influenced by Levy noises is due to various facts. Firstly, there is 
an intrinsic theoretical significance in the generalizations of previous ratchet models to 
account for cases where fluctuations present long tailed probability distribution, giving 
rise to anomalously large particle displacements. Moreover, since Levy noises may in- 
duce divergencies in the moments of the velocity distributions [5] , there is an additional 
challenge in providing the appropriate quantities to measure currents [SI EH] , particle 
dispersion [5[ [11] and efficiencies [9] . From another point of view, Levy ratchets may be 
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of direct interest for problems in magnetically confined fusion plasmas [5J HJ Q3] . Levy 
ratchets might be also of interest for atomic transport in dissipative optical lattices. 
Cold atom ratchets have been indeed recently studied experimentally [I5j [161 E3 HEl HE] 
and theoretically [201 EI], whereas anomalous dynamics described by Levy statistics 
were also reported for the same system [221 123] . 

In this paper we study the dynamics of a ratchet system influenced by two different 
signals inducing non equilibrium conditions: a dichotomic multiplicative noise and an 
additive Levy forcing. We analyze both the Langevin and Fokker-Plank approaches, 
focusing mainly on the latter. In different limit situations, the model generalizes various 
previous systems found in the literature. As we will see, for fast, slow and null flashing 
we recover different versions of the static Levy ratchet, while for a — > 2 we recover 
a version of the well known flashing ratchet with thermal noise [HE]- In particular 
we find that flashing allows both to enhance and diminish in a broad range the static 
ratchet current depending on the frequencies of flashing. Our results thus extend those 
for dichotomous flashing ratchets with gaussian noise to the case of broadly distributed 
noise. 

It is worth mentioning that a dichotomous flashing ratchet with Levy noise was 
shortly analyzed in the past as part of other studies on the influence of supperdiffusion 
on directional motion [21]. However, only heuristic arguments were given by considering 
an adiabatic approximation which assumes relatively slow flashing. Interestingly, as the 
study is previous to the findings in [5] and [6] , the contribution to directional transport 
during the "on" stage of the flashing was naturally ignored. 

The organization of the paper is as follows. In section [2] we introduce the model. 
In section [3] we give the solution of the Fokker-Planck equation. In section H] we analyze 
some limit situations which connect our results with previous results for other systems 
studied in the literature. In section [5] we present a detailed numerical analysis of the 
results considering a particular but standard ratchet potential. Section [6] is devoted to 
our conclusions and some final remarks. 



2. Model 



We consider the one-dimensional overdamped motion of a particle at position X(t) 
subject to a randomly fluctuating ratchet potential f(t)V(x), and to an additive Levy 
random force £(t). This situation is described by the Langevin equation 
d X 

±£ = -f(t)v'(x) + z{t). (l) 

We consider £(£) as a symmetric a-stable Levy noise with 1 < a < 2 and intensity x, 
described by the characteristic function (exp[ik J t * +A * £(t')dt']) = exp[—x\k\ a At] [SlEllT]. 
We consider a ratchet potential with spatial period L and amplitude AV, and the 
multiplicative noise f(t) as a dichotomous Markov process [251 I2S] taking the value 
of 1 (mode A) or (mode B), switching between these modes at rates Tab and Tba, 
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respectively. Equation flTj), thus, models a particle in a ratchet potential V(x) that 
switches between 'c nd T BA respectively 

(see Figure [[]) . We 'ansport that can be 




Figure 1. Scheme of the Levy ratchet model with dichotomic random flashing. 

generated in this system. This can be quantified by the long-time limit of the average 
velocity (dX/dt), where (...) stands for the average over all possible realizations of 
and /(£). Note that we only consider the domain 1 < a < 2 since it was shown that 
transport is highly inefficient in Levy ratchets for a < 1 due to the divergencies of the 
mean value of Levy distributions [SJE]. Moreover, the analysis in such parameter region 
demands special care for defining the mean velocity and current [5]. Note that, our 
parameter region of interest includes the case a = 2, for which the Levy noise becomes 
Gaussian white noise and (0Q) describes the dichotomous flashing ratchet studied in [27J. 
On the other hand, for T AB = 0, equation (JTJ corresponds to the static or non flashing 
Levy ratchet mentioned in the introduction and recently analyzed in [5j E]. 

The Langevin equation ([1]) was solved by using a standard algorithm (see [5] and 
references in j5] and [6]) for which the numbers distributed according to Levy laws were 
obtained by the method explained in [28J. 

As mentioned in the introduction, most of our analysis will be done using the 
Fokker-Planck approach to the problem. In order to present it, we first consider a 
more general situation in which the potential in state B can be non vanishing, and we 
introduce the probability densities P A = P A (x,t) and P B = P B (x,t) of finding the 
particle at position x and time t in the A and B modes respectively. For such a system 
we have the following coupled Fokker-Planck equations for the probability densities 

d t P A = d x (P A d x V A ) + X d%P A + T BA P B -T AB P A 

d t p B = d x (P B d x v B ) + xd«P B + r AB p A -r BA p B . u 

Here, d° = <9 a /<9|x|° stands for the Riesz-Feller fractional derivative of order a |29j, 
while V A (x) and V B (x) are the external potentials considered in states A and B. 
For V (x) = 0, the system ([2]) is equivalent to the Langevin equation ([T|). We will 
limit to such a case when performing our analysis of transport and presenting our 
numerical results. However, the method of solution of the Fokker-Planck equation will 
be introduced (up to some point) for an arbitrary potential V B (x), since this leads to 
more general and symmetrical expressions. 

Within the Fokker-Planck framework, the probability current J(x, t) = 
((dX(t)/dt)5(x — X(t))) is obtained from the continuity equation d t P = —d x J, where 
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P = P A + P B = (S(x — X(t))) is the probability density regardless of the mode. In 
the long-time limit d t P = d t P A = d t P B = and the current reaches the spatially 
constant value J characterizing the average directed motion or ratchet response. Since 
the steady-state solutions we are interested in, P A (x) and P B (x), are periodic with the 
period L, we can direclty solve (121) for dtP A = dtP B = in the interval < x < L with 
periodic boundary conditions. The steady-state current is then simply related to the 
average velocity in the Langevin description, J = (dX(t)/dt) /L. 



3. Solution of the Fokker-Planck Equation 

We exploit the fact that the steady-state solution we seek is periodic with the period 
L of the ratchet potential by using discrete Fourier transforms, f(x) = J2 q fg exp(igx), 
where q = 2Trn q /L with n q an integer and f k = J^(dx/L) f(x) exp(— ikx) denotes the 

correponding Fourier amplitude. In particular, the prescription [d"f(x)] k = — \k\ a fk for 
the Fourier transform of the fractional derivative, valid in an infinite or periodic support, 
greatly simplifies our method as it avoids the complications due to the nonlocal nature 
of the fractional Laplacian operator that can arise for different boundary conditions that 
break translational invariance in bounded domains [30]. Fourier transforming ([2]) gives 



dtP k A = -kZ q qV q A P A _ q -x\k\ a P k A + r BA p k B -T AB p k A 

dtP k B = -kZ q qV B P k \- X \k\ a P k B + T AB P A -T BA P B , 



where we have used (and use from now on) indistintevily the wave vector name, eg k, 
to denote both its value, k = 2irn k /L and its associated integer, n k , to avoid excesive 
notation. By construction equation ((3]) only admit periodic initial conditions with period 
L. This is not important as we will be interested in the steady-state limit d t P k ,B = 0, 
whose solution is unique and thus independent of the initial condition. Equation ([3]) 
must be solved with the constraint P A + P B = 1, arising from the normalization of 
the entire probability function J Q dx P = 1. This leads inmediately to the zero-mode 
steady-state solution 

Pf = ^ba/ (Xba + Tab) ^ 
Pq = Tab/ (Tba + Tab), 

which physically simply states that the probability of finding the system in the A and 
B modes regardless of the particle position, are controlled by the average fraction 
of time spent on each mode, t a /{t a + r B ) and t b /{t a + t b ) respectively, where 
ta/t b = T BA /T AB . 

Using the zero-mode solution of (jlj) the system ([3]) can be written as an 
inhomogeneous system for all the components P k and P k with k ^ 0, 

-kJ2^ k qV A P A _ q - X \k\ a P k A + T BA P B -T AB P A = k% A 6 
-kE q ^QV B P B _ q -x\k\ a P k B + T AB P A -T BA P B = eV B (l-6), {) 

where we have defined the duty ratio or fraction of time spent in the A mode, 
5 = T BA /(T BA + T AB ). In the particular case of a potential that is stochastically 
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switched on and off, we have V (x) = V(x) and V B (x) = and thus V B = 0. Using 
the second equation in the system of Eqs. (j5]) we obtain, 

pB = TaB pA (6) 

k x\k\ a + T BA k K J 

and therefore the problem is reduced to solve, 

- k '£qV q pA q -x\k\ a (l+ .J™ )p k A = ev k 5 (MO) (7) 
q+k \ X\k\ +LbaJ 

with P^ = 5 = 1- P,f , from (0}. In a more compact matrix notation we have 

(M + B)P A = CS, (8) 
where the components of the matrices M, D and the vector C are defined as, for k ^ 0, 
M kq = -k{k - q)V k _ q , 

D ^ = -x\k\ a (l + ^^)s kq , (9) 
C k = -M k0 = k 2 V k . 

In this way, M and C depend only on the Fourier components of the ratchet potencial 
V(x), while the diagonal matrix M concentrates, for fixed 5, all the dependence with 
the external parameters Tab, T b a, X an d ot- 

Finally, in order to compute the current we consider the Fourier-transformed 
continuity equation 

d t P k (t) = -ikJ k (t), (10) 
where P k {t) = P k {t) + P k {t). Adding the two equations (jBJ) together, we can identify: 
Jk = 1^ P k- q ~ ix\k\ a - l sg<k){P k + P k )- (11) 

As indicated in the previous section, the current in the stationary regime is a constant. 
This can be seen from (1101) . which for dtP k (t) = implies J k = for k ^ 0. Thus, 
assuming a > 1, the stationary current is simply 

J=J = -iJ2qV q P A q . (12) 
q 

Note that the solution for the general case V B (x) ^ demands dealing with the 
whole system ([5]). Clearly, although the procedure is a little bit more intricate, the 
equations can also be written in matrix form and solved for the variables P k and P k , 
enabling the calculation of the current. 

4. Dimensional analysis, relevant parameters and limit situations 



Considering a potential V (x) of amplitude AV and period L, a straight forward 
dimensional analysis of (171) shows that the Fokker-Planck solutions for the probability 
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distribution and current are of the form 

fx T AB L 2 X L 2 - 
P = f [ T ;6,a, 



L AV AV 

J= ^ M '^^'^J— ' (13) 

where / and g are dimensionless functions. Thus, there are essentially four relevant 
parameters. Note that the scaled 'on' rate r baL 2 / AV can be considered as a relevant 
parameter instead of 5 or instead of r abL 2 / 'AV . In the following we will speak 
alternatively of the three parameters depending on the feature we want to stress. Note 
that for 5 = 1, (JTj) reduces to the equation for the standard, non flashing, Levy ratchet 

mm- 

The relaxation time for a particle in the potential V(x) is expected to be a relevant 
characteristic time of the system. According to the deterministic part of (pQ), considering 
a typical length L and a velocity (AV)/L, we get L 2 /AV as the typical relaxation time. 
It is interesting to note thus, that the scaled parameters r abL 2 / AV and Y B aL 2 / AV are 
simply measures of the switching rates in the time scale of the inverse of the relaxation 
time. Equivalently, L a /x can be interpreted as a super- diffusion time which in the 
solution of the Fokker-Planck equation appears weighted to the relaxation time. Note 
that in the Gaussian limit we get the usual expression for the diffusion time = L 2 /x- 

By modifying Y A b and T BA but keeping constant their ratio (i.e. keeping 5 
constant), we can analytically study the limits of infinitely slow and infinitely fast 
switching without changing the relative residence times. Using (JTJ) the calculation of 
these limits is very simple and instructive. 

4-1. Infinitely slow switching 

In the infinitely slow switching case we have Tab and Tba — > 0, keeping 5 constant. 
Replacing this in ([7]) and using (BJ, we can write 

-* qV q (P k A _ q /5) - x\k\ a (P k A /s) = k 2 v k (k + o) 

(Pf/6) = 1. (i4j 

These equations are identical to the non- flashing Levy ratchet equations for Pj^/S. If 
we call P st (V) the real-space static Levy ratchet solution in a potential V, we have, in 
this slow switching limit, 

j slow = Jstiy) 5, 

where the second equation follows from (Tl2l) . Intuitively, the idea is that the rates are 
so large that the system reaches the steady state when the potential is 'on' and when 
it is 'off', and having the transients between these modes a negligible contribution, 
the total current is just an average of the two corresponding steady state currents. 
According to ([7]) the result is expected to be a very good approximation when the 
following relations hold between the scaled parameters: Y abL 2 / AV xL 2 ~ a / AV and 
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T BA L 2 /AV < xL 2 ~ a /AV (i.e. V AB < x/£ Q and T Bj4 < We will check this in 



section [5] when analyzing transport for the case of a standard ratchet potential. 
4-2. Infinitely fast switching 

In the infinitely fast switching case we have Tab — > 00 and T B a —> 00, keeping S 
constant. Replacing this in ([7]) the bracketed term becomes 1/5. We can thus write, 



but in the renormalized potential V 5. Thus, the solution in this fast-switching limit is 



In this case, the current for the flashing Levy ratchet is the same as for a system where 
the same potential is always 'on'but with an intensity rescaled by 5. This is intuitively 
evident: the potential is being switched so fastly that the particle can only 'feel' the 
temporal average of the potential (i.e. V(x)5). According to ([7]) the limit results in 
( IT71) are expected to be valid for Tab 3> x/ L a and Tba 3> x/L a - We have verified this 
using a standard different iable ratchet potential, as we will show in section |5j 

Note however that, at variance with what happens in the slow limit, the fast limit is 
not valid for all the Fourier modes. The denominator of the bracketed term in ([7]) shows 
that the fast approximation is only valid for Fourier modes that satisfy \k\ a <C Tab/x- 
Nevertheless, if the potential is an analytic function the probability density function, 
P(x) should also be an analytic function, which implies that the coefficients of its 
Fourier expansion decay exponentially with k. Therefore we expect that the higher 
Fourier modes that are not well approximated in the fast limit do not change the overall 
limit of the probability density. 

4-3. Perturbative analysis for fast and slow switching 

Using a simple perturbation analysis in ([7]) we can obtain the lowest order corrections 
to the infinitely slow and fast switching limits described above, with Tab and T^ B the 
small parameters, respectively, for a fixed value of S. 

Approaching the infinitely slow limit a natural ansatz is to propose, P A (x) ~ 
Psiow( x ) + ^abW(x), with P^ ow given by ( IT5|) . Inserting this ansatz in (J7|) we get the 
following result for the fourier-transformed first order correction for k 7^ 0, 



For various particular potentials we confirm numerically that (fl8|) has solutions with 
finite currents. As we will see, this is also consistent with the slow switching behaviour 
found in the full numerical solution discussed in the next section. We therefore conclude 




Pfast=Pst(VS)5 
J fast = Jst(VS) 



(17) 



kJ2^ q W k - q - X \k\ a W k = [P s fj fe . 



(18) 



q^k 
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that in the absence of particular symmetries, the first order correction generically does 
not vanish. We thus predict 

J — J slow ~ r \b (19) 

for small enough Tab (Xba) at fixed 5, with J s i ow given by ( !T5|) . 

Similarly, approaching the infinitely fast switching limit a natural anzats is to 
propose P A (x) ~ P A ast (x) + T^ B W (x) , with Pf ast (x) given by (1TTI) . Inserting this ansatz 
in ([7]) we get the following result for the fourier-transformed first order correction for 

-kJ2<lV q W k _ g - X \k\ a (W k /5) = -[(5- 1 - l)x\k\ a ] 2 [Pf ast } k . (20) 

q+k 

As before, in the absence of particular symmetries this equation yields generically 
solutions with a finite current. This is also confirmed by direct numerical evaluation of 
the last equation for particular potentials, and is also consistent with the fast switching 
behaviour of the full numerical solution presented in the next sections. We thus predict 
a non- vanishing first order correction 

J ~ J fast ~ r A £ (21) 

for large enough Tab (Xba) at fixed 5, with Jf ast given by f|T7|) . 

It is important to note that the asymptotic behavior of the currents that we have 
found (Equations ( IT9|) and (I2T!) ) are also valid for the case of Gaussian noise (a = 2). 
However, they do not coincide with the corresponding limits found for a randomly 
flashed triangular ratchet [27] and for a periodically pulsated ratchet with a frequency 
Q [2] . In this last case, the corrections analytically found in the slow and fast limits are 
J ~ Q 2 and J ~ — Q~ 2 respectively. The discrepancy is most likely due simply to the 
difference between a stochastic and a deterministic switching of the potential. 

In the case of [27] , however, the discrepancy is more significative because the system 
analyzed is also a (dichotomous) stochastically flashed ratchet. Using an heuristic 
argument, Astumian and Bier provide an expression for the current as a function of 
what we have called Tab (using 5 = 0.5), which in the slow and fast limits behaves as 
J ~ T 3 / B and J ~ exp(— Tab)- One of the possible reasons for this difference could be 
that that, in the case of a non analytic potential, the Fourier coefficients do not decay 
fast enough for our arguments to be valid. Another possible reason could be that the 
qualitative argument in [27J is not correct in the limits of fast and slow switching. 



5. Results for a standard ratchet potential 



Now we consider a standard ratchet potential [T] 



ZTT 



sin (2ttx) + - sin (Attx) 



(22) 



with period L = 1 and barriers of amplitude AV = V max — V m i n ~ 0.35. For this choice 
of the potential, we analyze the dependence of the current on the parameters x and a 
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characterizing the Levy noise, and on the transitions rates Tab and Tba- The results 
of this section correspond, thus, to the system (JTJ considering the potential V(x) given 
in f l22l) or, equivalently, to the Fokker-Planck equation (j2J) with V A (x) = V(x) and 
V B (x) = 0. The results for the current given throughout this section correspond mostly 
to the Fokker-Planck formalism. However, in order to show the complete agreement 
between the two formalisms we also include Langevin results for some specific cases. 

Concerning the solution of the Fokker-Planck equation, the Fourier transform of 
V(x) is simply 



k 4vr 



Sk,k+1 + ^k,k-l + ^(<5fc,fc+2 



(23) 



where 6ij is the Kroneker Delta. Thus, the matrix M given in (Q is pent a- diagonal. 
Despite this simplification, the system (jSJ) can only be solved numerically, for what we 
typically consider a number of Fourier modes ranging from N = 1000 for a ~ 2 to 
N = 2000 for a close to 1. This guarantees that our solutions closely approximate 
those of the continum limit if the Levy noise intensity satisfy \ > T BA/\kmax\ a , with 
kmax — ZttN/L the maximum wave vector. 



5.1. Symmetrical transition rates 

We first consider the case of equal transitions rates between A and B states. This means 
that the potential is 'on'and 'off' the same fraction of time, on average. We thus have 
5 = 0.5 and we define T = Tab = 

In Figure |2] we show the current as a function of the noise power \ for several values 
of the transition rate T and two different values of a. As usual in most ratchet systems, 
in all the cases studied we see that the current attains a maximum at an optimal noise 
intensity, while it decreases to zero for low and large intensities. The optimal value of 
X depends only slightly on T. As we can see, it decreases in a factor of order 1/2 when 
changing T in six orders of magnitude. We can also see that it decreases with increasing 
a. In all the cases studied the optimal value of x remains in the interval between AV/10 
and AV. 

Concerning the weak noise limit x — >■ 0, our numerical results seem to indicate that 
the current behaves as J ~ Xi m agreement with the findings for non flashing Levy 
ratchets in [10J. However, we have not been able to obtain the exact analytical limit 
law within our Fokker-Planck approach. 

Note that, for the two values of a analyzed in Figure El we find that T = 2 produces 
the largest value of J for almost all values of x, while larger and smaller transition rates 
lead almost always to smaller currents. The results in Figure |3] confirm the existence of 
an optimum value of T maximizing J, which is almost independent of a and grows with 
X- 

Figure |3] also shows the validity of the asymptotic formulas found in section H] for 
the limits of small and large T, and of the perturbative analysis close to them. Another 
fact that can be observed in Figure |3] is that, for a fixed flashing mechanism (fixed T) 
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Figure 2. Current J as a function of Levy noise intensity for different transition rates 
(r = 2 x 1CT 3 , r = 2x 1CT 1 , T = 2 x 10°, T = 2 x 10\ T = 2 x 10 3 ). (a) stability 
index a = 1.4. (6) stability index a — 1.8. All the continuous curves correspond to 
Fokkcr-Planck results, while the open circles were obtained from Langevin simulations 
for Tab = 2 x 10°. 



and a fixed noise intensity, Levy noise leads almost always to larger values of J than 
Gaussian noise. 

Finally, in Figure H] we analyze further the dependence of current with the noise 
power and the stability index. In order to identify a relevant noise region, we plot J as 
a function of a and x f° r an intermediate (near optimal) value of the transition rate. 
For a small, the range of noise intensities that maximizes the current is centered around 
X ~ 0.4. As the stability index increases, the optimal noise intensity shifts to smaller 
values. At the same time, the optimal noise range narrows down as we move to a = 2. 
It is interesting to note that a given current value can be obtained by combining x an d 
a in different ways. For instance, it is sometimes possible that a small noise power with 
a large value of a gives the same current as a large noise power with a small a. 

5.2. Non- Symmetrical transition rates 

Now we study the general case of different transitions rates between 'on'and 'off' states. 

In Figure [5] we show the current as a function of 5 for different values of a and 
considering a relatively small noise intensity x = 0.05. We see that the current vanishes 
for 8 ~ 0. This is because, at that limit, the potential remains most of the time 'off'. In 
contrast, for 5 — > 1, the potential is 'on'most of the time and we get the static or non 
flashing Levy ratchet limit. In such situation, the current depends on a and vanishes 
only in the case a — > 2, when the system approaches the equilibrium situation of a non 
flashing ratchet with Gaussian noise. 

The different panels of Figure [5] show us that the transition probability Tab plays a 
significant role in the current behavior (i.e. 5 alone does not determine the dynamics). 
For intermediate values of Tab (Figures \Mb) and[5]^c)) the current has a maximum at an 
intermediate duty ratio for almost all values of a. In contrast, in the case of very large 
or very small Tab (Figures E(a) and [5]^)), except for a equal or very close to a = 2, the 
current is a monotonic function of the duty ratio, and reaches the maximum for 5=1. 
Note that the linear dependence of J on 5 observed in figure Elja) corresponds to the 
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Figure 3. Current as a function of the switching rate for different values of a 
considering x — 0.05 (a),x — 0.2 (b) and \ — 0-6 ( c )- In panels (a), (b) and 
(c), the symbols at the extreme values of V indicate the infinitely slow (Eq. [15]) 
and infinitely fast (Eq. I17j) limit approximations showing complete agreement with 
the exact solutions. Panel(d) shows (J — J s iow) vs. T in the slow switching range 
confirming the validity of the asymptotic formula of Eq. (|f 9j) . The segment indicates 
a linear dependence on T for reference. Analogously, the results for (J — J fast) in 
panel (e) shows the validity of the perturbative analysis at fast switching fEql2Tj). The 
segment indicates a 1/T dependence. Results in panels (d) and (e) are for \ = 0.05. 

limit of slow switching indicated in section HI Namely, J~i5x J st (V). 

Figure |6] analyzes the dependence of J on 5 considering a larger noise intensity, and 
sweeping wide ranges of a and Tab- We see that in most cases the maximum current 
is achieved close to 5 = 1. The exceptions occur for values of a close enough to 2 (for 
instance that on Figure M^c)) and intermediate values of T^s- 

We can summarize the results of our analysis of J as a function of 5 at fixed 
Tab as follows. When considering small enough \ or large enough a, the inclusion 
of an appropriate flashing mechanism improves the performance of the static ratchet. 
In contrast, for large x or small a, the largest value of J is obtained considering a 
slow flashing mechanism close the static ratchet (i.e. 5 = 1). We thus see that the 
currents depend separately and richly on both, Tab and 5, allowing us to both enhance 
or decrease the static Levy ratchet currents for fixed noise parameters. 
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1.4 1.6 
a 

Figure 4. J vs a and x for a fixed value of transition rate r = 2 x 10 _1 . 
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Figure 5. Current versus duty ratio for different values of a and fixed x = 5 x f 2 
considering T AB = 2 x fCT 3 (a), Tab = 2 x KT 1 (b), T AB = 2 x 10° (c) and 
r\4 B = 2 x 10 2 (d) 
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Figure 6. Current versus duty ratio for a fixed value of the noise power \ = 0-2 and 
different values of Tab- Calculations for a = 1.2 (a), a = 1.6 (&), and a = 1.8 (c). 

6. Conclusions and final remarks 

We have studied the combined action of a Levy additive noise and a random dichotomic 
flashing in a ratchet system. Our results provide a complete generalization of previous 
studies on 'non-flashed' Levy ratchets and on standard flashing ratchets with Gaussian 
noises. 

We have presented a complete analysis of the two- variable fractional Fokker-Planck 
equation associated to the system. Our Fourier treatment allowed us to convert 
analytically the system of partial differential equations in an infinite-matrix linear 
system that can be easily solved numerically considering an appropriate truncation. 
Moreover, we were able to provide analytical asymptotic laws for slow and fast flashing 
that behave respectively as J — J s i ow ~ T and J — Jf as t ~ 1/r, where T is the flashing 
frequency The solution of the Fokker-Planck equation was given for an arbitrary 
periodic potential, and indicated also for the case in which the system switches randomly 
between two different potentials. 

Considering a standard ratchet potential we have systematically investigated the 
behavior of the current as a function of the stability index of the Levy noise, the noise 
intensity and the 'on' and 'off' rates of the flashing mechanisms. The Fokker-Planck 
results for the current were checked by means of Langevin calculations. We have found 
that random dichotomic flashing can produce a rich behaviour of the ratchet current. 
It allows both to enhance and diminish appreciably the static Levy ratchet current 
depending on the magnitude and relative magnitude of the flashing frequencies, and on 
the Levy noise parameters. A general statement to remark is that for small enough 
noise intensity or large enough stability index, a flashing mechanism can enhance the 
current of the static ratchet. Another relevant result indicates that, for a fixed flashing 
mechanism, the Levy noise gives larger current than the Gaussian noise in almost any 
situation. 

Our work thus contribute with quite general results and procedures to the 
understanding of the transport mechanisms on ratchets. In particular, to the rapid- 
growing new field of Levy ratchets. 
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